Microwave response of a two-dimensional electron stripe 
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Electromagnetic response of a finite- width two-dimensional electron stripe is theoretically studied. 
It is shown that retardation and radiative effects substantially modify the absorption spectrum of the 
system at microwave frequencies, leading to a non-trivial zigzag behavior of the magnetoplasmon- 
polariton modes in magnetic fields, similar to that recently observed by Kukushkin et al. [Phys. 
Rev. Lett. 90, 156801 (2003)]. 
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I. INTRODUCTION 

Plasma oscillations in low-dimensional electron sys- 
tems (ES), - inversion layers, systems on the sur- 
face of liquid helium, quantum wells, wires and dots, 
were experimentally studied since 1976 1,2,3 and 
seemed to be well understood (see e.gi 4 : 5 : 6 : 7 : 8 : 9 : 10 ! 11 ). 
Far-infrared transmission, as well as Raman spec- 
troscopy experiments quantitatively confirmed theoret- 
ically predicted 1 ^ 1 ^ plasmon and magnetoplasmon spec- 
tra, both in two-dimensional (2D) electron layers, and 
in systems with lower dimensionality (wires and dots). 
However, recent experimental studies 1 ^* 1 ^ of the mi- 
crowave response of macroscopic 2D electron disks with 
the diameter ~ 1 mm revealed very interesting and fully 
unexpected features of the magnetoplasmon spectra. Re- 
sponse of relatively small disks (with the diameter < 0.1 
mm) was very similar to that of microscopic dots. The 
spectrum consisted of two sharp resonances correspond- 
ing to the bulk and edge magnetoplasmon modes with 

the frequencies u)±(B) = a/w 2 + w 2 /4 ± u) c /2, where w c 

is the cyclotron frequency and u p the plasma frequency at 
zero magnetic field B = 0. In larger samples, however, 
the upper (bulk-magnetoplasmon) mode exhibited sub- 
stantial changes. At a certain finite value of B this mode 
intersected the cyclotron resonance (CR) line u> = u> c and 
disappeared experiencing a dramatic broadening. A sec- 
ond, and in some cases a third mode, with comparable or 
even larger amplitude, appeared in the absorption spec- 
trum at higher frequencies and showed similar behavior. 
As a result, the -B-field dependence of the bulk magneto- 
plasma modes exhibited intriguing zigzag behavior as a 
function of magnetic field at u> > w p 14,15 . 

It was shown ini 4 * 1 ^ that the observed features of the 
magnetoplasmon spectra are due to the influence of re- 
tardation effects. At B = the measured resonance fre- 
quencies turned out to be in a good quantitative agree- 
ment with Refi 1 ^, and at finite B calculations did reveal 
a certain zigzag behavior of the resonant bulk magneto- 
plasmon modes. However the physical reasons of such a 
strange behavior of the modes at finite B, as well as the 
overall behavior of the linewidth of resonances, have not 
been explained. 

The goal of this paper is to theoretically study elec- 



tromagnetic response of a finite-size macroscopic 2DES. 
We consider a 2D stripe of a finite width W and 
focus on the regime, where the retardation and ra- 
diative effects become important (for mm-size 2D 
samples this corresponds to microwave frequencies - 
the parameter range where very interesting microwave 



experiments 



17,18.19.20,21.22,23.24 



have been recently done). 



Although the stripe geometry differs from the experi- 
mental one (disks), our calculations reproduce all the 
experimental findingsi 4 * 1 ^ and allow one to understand 
the physics of the observed effects. It will be shown that 
the unusual zigzag behavior of the bulk-magnetoplasmon 
modes (more exactly, magnetoplasmon-polaritons) is ex- 
plained by an interplay of radiative and dissipative con- 
tributions to their linewidth. 

In an infinite sample the spectrum of 2D magneto- 
plasmons has been considered (taking into account re- 
tardation but ignoring the damping of modes) by Chiu 
and QuinnH in 1974. Their results are applicable for 
the description of far-infrared response, when the size 
of the sample is large as compared to the magnetoplas- 
mon and radiation wavelength, and the issue of linewidth 
is not very important. We consider a finite-size ge- 
ometry, more relevant in view of the recent microwave 
experimenta 14 i 15 i 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 , where the radiation 
and the plasmon wavelengths exceed or comparable with 
the sample dimensions, and where the linewidth of modes 
essentially determines the observed non-trivial features of 
the absorption spectra. 



II. THEORY 

Assume that the 2D electron- gas stripe (— W/2 < x < 
W/2, — oo < y < oo) lies in the plane z = 0, the external 
magnetic field B = (0,0, B) is directed along the z-axis, 
and the dielectric constant e of surrounding medium is 
uniform in all the space. Let an external plane elec- 
tromagnetic wave with the frequency u> and the elec- 
tric field E ext (r,t) = E ext e*"(*^/ c - t ) be normally inci- 
dent upon the structure. It induces in the 2D system 
oscillating charge, current j, and - around the 2DES 
- electric and magnetic fields. The total electric field 
E tot = W xt + E md , given by the sum of the external 
and the scattered (induced) fields, satisfies the Maxwell 
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equation 

V x (V x E tot ) + - 2 -^~ = -^ Ft i(x,t)S(z). (1) 

Looking for a solution, proportional to e~ luJt , and as- 
suming that fields and currents do not depend on the 
y-coordinate, we rewrite Eq. JQ) in an integral form for 
the x- and y-components of the induced electric field 

■ ,00 f W/2 

E™ d (x) = / dq x Ke^ x , ' ./,■:./:,//. 

UJt J-00 J-W/2 

(2) 

/oo 7 pW/2 
^ e ^ x / e~ lq * x ' Jy (x')dx' . 
-00 K J-W/2 

(3) 



Here k — y^ — Lu 2 e/c 2 determines the z-dependence of 
fields (cx e _K ' 2 '). Dependent on the interval of the dis- 
integration, Eqs. ijSl - © describe both the evanescent 
induced fields (at q 2 > uj 2 e/c 2 ) and the scattered outgo- 
ing wave (at q 2 < u> 2 e/c 2 ); the sign of Im k should be 
properly chosen to satisfy the scattering boundary con- 
ditions at \z\ — > 00. 

Now we need to relate the current and the total elec- 
tric field inside the stripe. We do not consider in this pa- 
per possible nonlocal effects and use the relation j a (r) = 
Vaf3E t p t {r)6{W/2 — \x\), where a a/ 3 is the conductivity 
tensor of the 2DES, dependent on the frequency, mag- 
netic field and the momentum relaxation rate, and 9(x) 
is the step function. Together with Eqs. J2J © this 
gives a system of two integral equations for the electric 
current inside the stripe \x\ < W/2, 



jx(x) = <J xa E ex 



W/2 

dq x ne lq * x / e~ lq * x ' j x {x')dx' + 

W/2 



j c-W/2 

^e^ x / e~ iq * x ' j y (x')dx' , (4) 

K J-W/2 



dq x Ke iq * x 



W/2 



W/2 



e~ iq " x j x {x')dx' + 



yy 



dq x 



W/2 



W/2 



jy(x')dx'. (5) 



These equations are solved by reduction to a matrix form. 
Expanding the currents in a complete set of functions, 



jx{x) 



1 ,," i cos 



2ltX 

~w 



1 

n+- 



(6) 



and following the standard procedure we get the following 
system of algebraic equations 



2(-iy 



71=0 



n=0 



V M {21) A +VM (22) B - 

/ j lvl mn -rt-ri \ / J ul mn D n — 



n=0 

where 



n=0 
1V1 mn — 

M< 12 > = 

mn 



?r(m + 1/2) 

2(-l) m 
7r(m + 1/2) 



a E ext 



(7) 



°ya rj a ; 



Stt 2 ^ (+) 

mn ^ Tiy ' J mn 



2iua xy W M 

9 ,J inn \ iL )i 



8ir 2 ia y: 



(8) 

(9) 
(10) 

(11) 



M (22) - 5 



mn u mn 



2iua yy W ( _) 

2 ■'mn V ji 7j 



and = u)Wy/e/2irc. The integrals are defined as 



(12) 



4±)(n) = / dQ(Q 2 - O 2 ) ±1 / 2 7m (0)7n(Q), 
where 



7n(Q) 



(-l) n+1 (n + 1/2) costtQ 

w Q 2 -(n+l/2) 2 



(13) 



(14) 



The part of the integrals in l|13(l from Q = — f2 to Q = 
gives imaginary contribution to the matrix elements in 
0-© and physically corresponds to the radiative decay. 

Having solved equations 0-© and having got the co- 
efficients A n , B n in terms of the external field E ext , we 
can calculate the current, Eq. ©, and the fields, Eqs. 

inside the stripe. Then, calculating the Joule 

heat Q = J W ^ 2 dx(j*~E + jE*)/4 and dividing it by the 

energy flow S z = c^l(\E ext \ 2 + \E^ xt \ 2 )/8TT, we get the 
absorption cross section Q/S z (dimensionality cm). All 
plots of the next Section show the frequency, magnetic 
field and other dependencies of the dimensionless "ab- 
sorption coefficient" 
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Q 



s z w 7 



(15) 
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defined as the absorption cross section Q/S z normalized 
by the width of the sample. 



III. RESULTS 
A. Preliminary notes 

The plasma frequency in a 2DES, calculated in the 
quasi-static approximation is given by the formula^ 



where 



2im s e 
< 



(16) 



where n s , e and m* are the density, the charge and 
the effective mass of 2D electrons, and q is the plas- 
mon wavevector. In a stripe of a finite width W the 
lowest quasistatic eigen-mode frequency is expected at 
q w ir/W. Therefore we choose for the frequency unit 
the quantity 



UJ 



2'K 2 n s e 2 
m*eW 



(17) 



for a typical macroscopic (mm-size) GaAs 2D sample its 
value lies in the tens-of-gigahertz range, 



^ = 43.3 x 
2tt 



/n s [10 n cm" 



eW[i 



[GHz 



(18) 



(for the effective mass we use m*/m ~ 0.067). Every- 
where below we measure the frequencies in units u/loq, 
magnetic fields in units w c /wo, and the momentum relax- 
ation rate 7 = ejm* \jl in units 7/wo (m is the mobility), 
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0.1 



/ eVK[mm] 
uT Q ~ /x[10 6 cm 2 /Vs] V n s [10 n cm" 



(19) 



We also introduce the retardation parameter a, defined 
as the ratio of the frequency of the quasi-static 2D plas- 
mon 1)16(1 with the wavector q — ir/W to the frequency 
of light TTc/^/eW with the same wavector, 



ujQyfiW 2n s e 2 W 



(20) 



in the considered case of the uniform dielectric constant 
the retardation parameter a does not depend on e. Its 
numerical value is negligibly small for microscopic 2DES 
(for instance for quantum wires and dots), but can be of 
order unity for mm-size 2D samples, 



a 



cm 2 ] x W^mm] 
The formula H20fl can be also presented as 

r 

a = — , 

UJQ 



(21) 



(22) 



r = 



27rn s e 2 
m*C\/e 



(23) 



is the radiative decay rate in an infinite sample (see a 
discussion ini&). 

We solve the system of equations {Zl~© numerically, 
restricting ourselves by a finite number Nh of harmonics 
in the expansions ©■ The matrix equation 0-0 then 
has the size 2 Nh x 2Nh ■ Convergence of the solution with 
respect to Nh has been checked. All results shown below 
have been obtained with Nh — 16. For the conductivity 
tensor a a p we have used the Drude model 
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' mi 



m* (lo + ij) 2 — uj 2 



n s e* 



777* (u) + ij) 2 — id 2 



(24) 



(25) 



B. Zero magnetic field 

We begin the presentation of results with the case of 
B = 0. We assume that the incident electromagnetic 
wave is linearly polarized, and show results only for the 
perpendicular polarization (Ey Xt = 0). The case of the 
parallel polarization is less interesting for the goals of this 
work. 

Figure exhibits the absorption spectra of a stripe as 
a function of frequency, at 7/^0 = 0.05 and a few val- 
ues of the retardation parameter a. One sees that at 
small a (the quasistatic limit) the absorption spectrum 
is similar to that of quantum wires: a strong peak corre- 
sponding to the fundamental 2D plasmon mode with q ~ 
ir/W is accompanied by a number of very weak higher- 
harmonics peaks, corresponding to q ~ (277 — 1)tt/W 
with n = 2, 3, . . .. Due to the symmetry of the exter- 
nal electric field E^. xt (x) = E%. xt = const, the even 2D 
plasmon modes with q ~ 2m: /W are not excited. When 
a increases, the fundamental mode first increases in am- 
plitude, reaches its maximum, and then decreases, ad- 
ditionally experiencing a very strong broadening. The 
behavior of the second mode is similar, but it reaches its 
maximum at a larger value of a. As a result, at a certain 
a the second-mode peak has a larger amplitude and a 
smaller linewidth than the first one (see e.g. the curve for 
a = 0.8). The third, fourth and all other modes behave 
similarly, but each subsequent mode reaches its maxi- 
mum (and becomes the dominant one in the spectrum) 
at a higher value of the retardation parameter. In gen- 
eral, as seen from the Figure ^ at a > 1 the absorption 
spectrum has the shape qualitatively different from that 
in the quasistatic limit: it consists of many peaks with 
quite comparable amplitudes, in contrast to the quantum 
wire/dot spectra, when the higher modes are hardly vis- 
ible and only the lowest one is dominant. All the modes, 
as seen from Fig. ^ experience quite essential red-shift, 
in agreement withiSiiii^. 
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FIG. 1: Absorption spectra of a 2D stripe at 7 /ljo ~ 0.05 and 
at different retardation parameters. 16 curves are shown in 
Fig. (a) with a varying from 0.1 up to 1.6 with the step 0.1. 
Fig. (b) provides a more detailed view on the same spectra 
around the first and the second modes. 



C. Finite magnetic fields 

Figure [3 illustrates the behavior of modes in finite 
magnetic fields. If a is small, Figure the absorp- 
tion spectrum is, again, very similar to that of quantum 
wires. The fundamental mode approaches the CR line ac- 
cording to the simple law uj mp (B) = 



] + LUc, slightly 
increasing in the linewidth and decreasing in the am- 
plitude. All the higher modes behave similarly; no one 
of them becomes stronger than the fundamental mode. 
Such a behavior was many times observed in experiments 
on quantum-wire arrays, see e.g&l&. 

If a becomes larger, Fig. the lowest, fundamen- 
tal mode losses its amplitude very quickly approaching 
the CR line, and the second mode becomes stronger. At 
even higher frequencies, and at larger values of a, Fig. 

the same happens with the second mode: it disap- 
pears, approaching and intersecting the CR, and releases 



its strength to the next mode. Experimentally, this man- 
ifests itself as a zigzag behavior of resonances^i^ (we re- 
mind that in experiments the absorption was measured 
as a function of magnetic field at a fixed microwave fre- 
quency): when lo increases, the lower- mode peak gets 
broader and disappears, intersecting the CR line, while 
a new, stronger peak, corresponding to the next, higher 
mode, arises at smaller B. At even larger values of ev, Fig. 
01, such a relay-race of oscillator strengths between the 
modes becomes more and more pronounced, especially 
for higher modes. The frequency of the lowest modes ex- 
hibits a very weak magnetic-field dependence at large a, 
Fig. Hi. 

It should be emphasized that the predicted features 
can be observed only in samples with a sufficiently high 
electron mobility. Figure [3] exhibits the absorption spec- 
tra at a = 0.4 (like in Figure EJj), but for a 4 times 
higher scattering rate (7/0*0 = 0.2). One sees that the 
non-trivial exchange by oscillator strengths between the 
different modes almost disappears and the whole spec- 
trum reminds now the one shown in Fig. (with broader 
linewidths). 



D. Discussion 

As known, in an infinite 2DES the CR has a maximum 
when the radiative and dissipative decay rates equal to 
each other. The absorption coefficient in this case (for 
a wave circularly polarized in the CR direction) has the 
forr 



,16,25 



.4. 



2 7 r 



(a* - o c )2 + ( 7 + Tf 



(26) 



which shows that the CR peak decreases and tends to 
zero both at T/ 7 — > and T/ 7 — > 00, and reaches its 
maximum (= 1/2) at r/7 = 1. Similar situation is the 
case in a finite-width stripe, and it is the competition 
between the radiative and non-radiative processes that 
explains the zigzag behavior of modes. 

Consider first the case B = 0, Fig. In the quasistatic 
limit a — > the linewidth of the n-th mode is determined 
by 7, while its amplitudes - by the radiative decay rate 
r„, which gradually decrease with n because of the grad- 
ual decrease of the corresponding dipole matrix elements 
(we stress that r ra are not identical to the radiative decay 
rate T (|23|l in an infinite sample). The growth of a leads 
to the increase of the radiative losses. When the radia- 
tive decay rate r„ of the n-th mode becomes equal to 
the dissipative decay rate 7, the amplitude of this mode 
reachs its maximum. As T n is a decreasing function of 
n, the absorption maxima for the higher modes are ob- 
served at larger values of a. This explains the behavior 
of resonances at B = 0, Fig. ^ 

In order to understand the behavior of modes at finite 
magnetic fields, it is helpful to compare our results with 
those of RefiH. Fi gure 0] exhibits the calculated ini£ 
dependence of frequency on magnetic field. As ini£ the 
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FIG. 2: Absorption spectra of a 2D stripe as a function of frequency and magnetic field, at "//ujo = 0.05 and at a few 
representative values of the retardation parameter a = 0.1, 0.4, 0.8 and 1.4. Each plot contains 21 curves for lo c /luo varying 
from to 4 with the step 0.2. 



eigen-value problem with a continuous magnetoplasmon 
wavevector q was considered, we normalize (only on this 
plot) the frequencies uj and u> c to aj p (q), Eq. I|16|) . an d 
use the retardation parameter defined as 

a^a%) = ^^ = 4 T ; (27) 
cq LU p (q) 

Lu p (q) and a*(q) are reduced to loq and a, if g = 7T/VF. 
As seen from Fig. the magnetoplasmon dispersion ex- 
hibits two different regimes. If Lo c juj v <c 1/cr*, it follows 

the quasistatic behavior u> — yk>% + wj? . The interaction 

of magnetoplasmons with light in this regime is small. If 
<-<Jc/^p ^ the magnetoplasmon mode intersects the 

CR line and tends to the asymptote oj/uj p — 1 /a*, which 
corresponds to the dispersion of light u> = cq/^fe. In 
this regime, the magnetoplasmon strongly interacts with 
light. In the infinite-sample geometryi^ the interaction 
of plasma modes with light thus manifests itself in the 



change of the dispersion of modes. 

In the finite-size geometry, strong interaction of mag- 
netoplasmons with light in this regime manifests itself 
in the strong radiative decay of modes. As in a stripe 
the n-th mode corresponds to q = n(2n — 1)/W, the 
strong interaction regime should be expected at lu/ujq ~ 
(2n — l)/a. As seen from Fig. |2 a substantial broaden- 
ing of the modes is observed, indeed, when the n-th mode 
frequency lo/loq approaches (2n — l)/a. As is also seen 
from these considerations, the modes with smaller n ap- 
proach the strong interaction regime and die at lower fre- 
quencies, releasing their oscillator strength to the higher 
modes. This explains the non-trivial zigzag behavior of 
the magnetoplasmon-polariton resonances, which is seen 
in Fig. [21 and was experimentally observed in Refsii^i^. 

We complete our discussion by answering the question, 
what happens with the spectra of Figure[21in the limit of 
a large sample. How the complicated structure of modes, 
shown in Fig. [21 is transformed to a simple, single-mode 
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FIG. 3: Absorption spectra of a 2D stripe as a function of 
frequency and magnetic field, at 7/wo = 0.2 and a = 0.4. 




-+= 1 1 1 1 1 1 1 1 1 1 1 1 

1 2 3 4 5 6 

co /co (q) 

FIG. 4: The frequency versus magnetic field dependence of 
magnetoplasmon-polariton modes, as obtained in Ref^. The 
units Lo/u} p (q), u] c /oj p (q), and a* — T/u p {q) are used on this 
plot; for explanations see the text. At large uj c /u) P (q) the 
curves tend to the asymptotes a>/u> p (q) = l/a*(q), which cor- 
respond to the dispersion of light uj = cq/yft. 

behavior of Eq. I|26l) in the limit W — > 00? To an- 
swer this question, notice that our frequency unit loq is 



proportional to 1/VW and decreases when W — > 00. 
This means that the dimensionless dissipation parameter 
"//luq grows up, and that the multiple-resonance region 
observed at u> ~ ujq is shifted to smaller frequencies in 
real (dimensional) units. As seen for example from Fig. 
in this limit one does get a single-resonance behavior 
of the absorption spectrum. 

IV. CONCLUSIONS 

As seen from our analysis, the most interesting effects, 
related to retardation, can be observed in 2DESs under 
the condition 

7 < lu < r. (28) 

This implies, first of all, that the parameter 7/r should 
be small. The condition 7/r <C 1 imposes the upper 
boundary on the zero-B resistance of 2D samples: it 
should be much smaller than the wave resistance of vac- 
uum, 377 Ohm (it is so, if the electron mobility in typical 
GaAs samples exceeds ~ 10 6 cm 2 /Vs). Second, Eq. i(2"5|) 
implies that the typical size of the samples should lie in a 
certain window, which extends (for typical material pa- 
rameters) from about 1 mm up to a few cm (dependent 
on the mobility). Recently published very interesting mi- 
crowave experiments on the 2DES 14 i 15 i 17 i 18 . 19 i 20 i 21 i 22 i 23 i 24 
have been done in this range of parameters, and hence 
the phenomena discussed here should be important for 
their proper interpretation 1 ^. 

To summarize, we have theoretically studied mi- 
crowave absorption spectra of a finite- width 2D stripe, 
taking into account retardation effects. We have shown 
that retardation substantially modifies the microwave re- 
sponse of macroscopic 2DESs, especially in the very-high- 
mobility samples. Our results are in very good agreement 
with experimental observations, and explain the myste- 
rious zigzag dispersion of magnetoplasmon modes, dis- 
covered in 1 * 1 ^. Further experimental studies of the pre- 
dicted features in a broader range of parameters would 
be very interesting. 
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